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Abstract 

We calculate the electric-dipole and magnetic-quadrupole form factors of the deuteron that arise 
as a low-energy manifestation of parity and time-reversal violation in quark-gluon interactions. 
We consider the QCD vacuum angle and the dimension-six operators that originate from physics 
beyond the Standard Model: the quark electric and chromo-electric dipole moments, and the gluon 
chromo-electric dipole moment. Within the framework of two- flavor chiral perturbation theory, we 
show that in combination with the nucleon electric dipole moment, the deuteron moments would 
allow an identification of the dominant source(s) of symmetry violation. 
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Permanent electric dipole moments (EDMs) of particles, nuclei, atoms, and molecules 
violate both parity (P) and time-reversal (T), or equivalently CP, invariance [IJ. Since CP 
violation due to quark mixing in the electroweak sector of the Standard Model (SM) seems 
insufficient to explain the observed matter-antimatter asymmetry in the universe p] and 
predicts immeasurably small values for the EDMs of nucleons and nuclei ji3|, searches for 
nonzero EDMs are an excellent probe for new sources of CP violation. Experiments with 
ultracold neutrons are in preparation that aim to improve the bound on the neutron EDM |1] 
by two orders of magnitude [3]. Moreover, plans exist to measure the EDMs of light ions, 
in particular the proton and the deuteron, in storage ring experiments at similar levels of 
accuracy The current generation of ongoing and planned EDM experiments probes the 
same energy scales as the LHC, and there are strong expectations that nonzero results will 
soon be found. 

An outstanding theoretical issue is to identify the fundamental CP-violating mechanisms, 
and in particular to relate a positive signal in the EDM experiments to the P- and T- 
violating {fl") sources at the quark-gluon level. The SM contains the fl" QCD 6 term 
[7], which has dimension four and would be expected to give the main contribution to 
hadronic fl". However, since the experimental upper limit on the neutron EDM constrains 
6 to be unnaturally small [8], ^< 10~^°, possible contributions from higher- dimensional fl" 
sources can be relevant, or even dominant. These higher-dimensional operators have their 
origin beyond the SM, in an ultraviolet complete theory at a high energy scale M^, for 
example a supersymmetric version of the SM [H]. The first such effective operators one 
encounters have dimension six [101 [H], viz. the quark EDM (qEDM) and the quark and 
gluon chromo-electric dipole moments (qCEDM and gCEDM). We show that in combination 
with the nucleon EDM, a measurement of the deuteron fl" form factors (FFs) would allow 
the disentanglement of these sources. 

The difficulty in calculating such low-energy observables stems from the breakdown of 
perturbation theory in the QCD coupling constant below the characteristic QCD scale 
^QCD ~ 27rP7r — 1-2 GeV, with F^, = 185 MeV the pion decay constant. For processes 
at momenta Q ~ rnT^, the mass of the lightest hadron, the pion, we can nevertheless ex- 
press observables in a controlled expansion in powers of Q using chiral perturbation theory 
(ChPT) [12] with the two lightest quark flavors u and d. This effective fleld theory (EFT) 
involving pions, nucleons and photons correctly incorporates the (approximate) symmetries 
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of QCD, in particular the spontaneously and explicitly broken 5*0(4) cliiral symmetry, and 
describes low-energy physics in a model-independent way. All effective hadronic interactions 
that transform under symmetries as terms in the QCD Lagrangian are allowed, each one 
being associated with a parameter, or low-energy constant (LEG), which can be estimated 
using naive dimensional analysis (NDA) [HI [13]. Since the sources break chiral symme- 
try in different ways [H], they can in principle be distinguished by the form and expected 
strength of their hadronic interactions. 

The nucleon EDFF partially reflects the source at the quark-gluon level [T3HT7]: a 
measurement of the nucleon EDM and its Schiff moment [ISj could distinguish between 
originating from the 9 term or qCEDM on the one hand, and the qEDM or gCEDM on the 
other. Since also the deuteron can be analyzed with firm theoretical tools, we focus here on 
its electromagnetic FFs which, part from small relativistic corrections, are defined from 
the part of the electromagnetic current, J^, by 

{p\3\Jh\P^^) = -e'''q'Fn{f), (1) 
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where \p, i) denotes a deuteron state of momentum p and polarization 6^ in the rest frame, 
normalized so that {p',j\p,i) = a/1 + p^/m^(27r)^5'^^)(g q = p — p' is the photon mo- 
mentum, K = {p' + p)/2, and = 2m n — + . . . is the deuteron mass in terms 
of the nucleon mass rriN and the binding momentum 7. We show that a combination of 
the deuteron EDM, = -Fd(O), and nucleon EDM can separate the qCEDM from the 
other fl" sources, and that a measurement of the magnetic quadrupole moment (MQM), 
■Md = -^m(O), is sensitive to the 9 term. 

P- and T-conserving (PT) pion interactions are relatively weak, because they proceed 
through derivatives when originating from the couplings of quarks and gluons, which are 
chiral-symmetric, or through powers of the small chiral-symmetry breaking parameters, the 
quark masses niu^d ^ rh = C(m^/MQCD) and the proton charge — e = ^/4^la^. In the one- 
nucleon sector the expansion is in powers of Q/Mqcd, but subtleties in the power counting 
arise in the two- nucleon (NN) sector [19] , where one has to accommodate the large scattering 
lengths in the ^5*0 and ^5*1 waves, and the related presence of the unnaturally shallow bound 
state in the ^Si wave, the deuteron. Such fine-tuning can be incorporated if the AW LECs 
are assigned a scaling [20] with inverse powers of the small scale 7 ~ 45 MeV. Moreover, the 



strength of pion exchange among nucleons is set by Mtvat ~ AnF^/mN — 450 MeV, and for 
momenta Q < Mat/v pions can be treated perturbatively [20] , the deuteron arising when the 
leading NN contact interaction is summed to all orders. NN observables are amenable to an 
additional expansion in powers of Q/Mnn, although the breakdown scale of this expansion 
in scattering |20] suggests that Mnn is smaller than guessed above by a factor of 2 or 3. 
This scheme was used to describe the low-energy properties of the deuteron, in particular 
its PT (charge and magnetic dipole) and fT (anapole) FFs pi] . We calculate here the 
FFs in leading order (LO) for the first time. 

At the level of the quark field q = {u dY and the photon and gluon field strengths F^^ 
and G^j^, the lowest-dimension sources are 

Cjfrp = m^e qi-f^q - -q {do + d^Ts) a^^-f^q F^^ 

-y {do + rfars) a^''l,X''q G^, + ^e^'^'^ T'^G^^^G'/GI^, (3) 

in terms of the Pauli isospin matrices r, the Dirac spin matrices 75 and a^'^, and the 
Gell-Mann color matrices A*^ and structure constants f"'^'^. The first term, with = 
fnufTT'd / {f^u + f^d) — fn/2, incorporates the QCD angle i9 [71 [8]. The second (third) term 
represents the isoscalar do {do) and isovector ^3 (^3) components of the qEDM (qCEDM). In 
the last term, dw is the gCEDM [H]. Because of electroweak SU (2) x U (1) gauge symmetry, 
the qEDM and qCEDM are proportional to the vacuum expectation value of the Higgs field, 
and therefore have effective dimension six \iQl, being suppressed by two powers of Mr. We 
write [16] 

*-of^V ''H-ofeV (4) 



where 5, 5, and w are dimensionless parameters that depend on the mechanisms of elec- 
troweak and PT breaking, and on the running from the electroweak scale Mw to low energies; 
their sizes can be calculated in specific high-energy models in terms of coupling constants 
and complex phases [HIH]- Our approach is limited to low energies, where the contributions 
associated with heavier quarks can be buried in the LECs. Effects of higher-dimension 
sources should be suppressed by M^^/M^,. 

While all interactions in Eq. ^ break P and T, each transforms under 5*0(4) in a 
characteristic way [TH4T7] . The 6 term is the fourth component of the same SO (4) vector P = 
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{qTq^qi'^^q) that leads to isospin breaking [Tll[T5l[T7], and thus generates EFT interactions 
that transform as fourth components of 5*0(4) vectors made out of hadronic fields, with 
coefficients related to those of PT interactions. Similarly, the qCEDM and qEDM both 
break chiral symmetry as combinations of fourth and third components of two other 5*0(4) 
vectors [161 lEj- The gCEDM does not break chiral symmetry. As in the PT case, we use 
NDA to estimate the strength of the effective interactions. 

The relevant Lagrangian [H] in terms of the pion field tt and the heavy-nucleon field 
N = {p n)^ oi velocity and spin is 

+|^e^'^°^t;„ NSp-K ■ tNF^^ + Oq [NN d^" [NS^N) - NtN ■ d^" {NtS^N)] 

f^SpN NSxN d^F^, + ..., (5) 



where e°^^^ = 1 



are pion-nucleon couplings, 



contributes to the short-range isoscalar nucleon EDM 
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'^^ = 0\e6^^^,\ (8) 



is a pion-nucleon-photon interaction. 



MqcdM^ 
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C, = 0\w 7"" (9) 



is the leading ff NN contact LEG, 
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parametrizes short-distance NN currents, and ". . ." stand for terms that only contribute 
to the FFs at higher orders. For i9, the link with isospin violation [H] implies — 5^Af(^d+ 





FIG. 1: LO diagrams for the deuteron EDFF (a, b) and MQFF (c, d, e, /). Solid, dashed, and 
wavy hnes represent nucleons, pions, and photons. A square marks f^, and the other vertices PT 
interactions: leading (filled circles) and subleading (circled circles). The hatched vertex represents 
the deuteron state. Only one topology per diagram is shown. 

mu)0 /2{md — mu) ~ 2.86 MeV, using lattice QCD input [22] for the quark-mass piece of 
the nucleon mass difference 6m n ■ 

The calculation of the EDFF and MQFF involves at LO the diagrams of Fig. [TJ where the 
squares denote interactions from C^. The circles denote well-known PT interactions, see 
e.g. Ref. [23]. The pion-nucleon vertex is the standard axial- vector coupling, = 1.27. The 
photon vertex denoted by a filled circle is the coupling to the charge e, and that denoted by 
a circled circle is the magnetic coupling parametrized by the anomalous magnetic moments, 
the isoscalar kq = —0.12 and the isovector ki = 3.71. The hatched circles denote deuteron 
states [21] obtained from the iteration of the leading NN contact interaction, whose LEG 
can be eliminated in favor of 7. We use dimensional regularization with power-divergence 
subtraction [20] at a renormalization scale /z. Our results depend on the ratio ^ = 
and on three functions of the momentum in the ratio x = |g I/47: 



which originates in a bubble with one photon coupling and appears also in the charge FF [21] , 
and two complicated functions that result from two-loop diagrams with a pion propagator, 
which can be expanded as 




(11) 




(12) 
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The scale of momentum variation is set by 47. 

The LO deuteron EDFF is due to diagrams a and 6, 

where the first term is dominant for 6, qEDM, and gCEDM, and the second one for qCEDM. 
The LO MQFF comes from diagrams c, d, e, and /, 
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go{l + Ko) + y(l + Ki] 



:i + 2{) 
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+ - ^) ^ + T:Ff (^rT^^3(x) + 21n , (15) 

where at this order (Jq originates from 6 and qCEDM, gi from qCEDM only, Cq from gCEDM, 
and M and c^r from qEDM. 

We can now discuss the implications of the various sources for the deuteron EDFF 
and MQFF. In Table |T] we list the orders of magnitude for the deuteron EDM, dd, the ratio of 
deuteron-to-neutron EDMs, dd/dn, and the ratio of the deuteron MQM and EDM, Aid/dd, 
for the different sources. Just as for dn [HI [151 [16], ^ dd signal by itself could be attributed 
to any source with a parameter of appropriate size. For 6, qEDM, and gCEDM the deuteron 
EDFF is determined by the LO isoscalar nucleon EDM, and thus well approximated by the 
sum of neutron and proton EDM. For 6 in particular, using the most important long- 
range contributions, which appear at NLO, as a lower bound for do [171 [21], one finds 
\dd\ ^ 2.8 ■ 10~^ 6 efm. If, however, the dominant/^ source is the qCEDM, dd comes mainly 
from neutral-pion exchange. A measurement of \dd\ significantly larger than \dn\ would be 
indicative of a qCEDM. A null-measurement at the 10~^^ efm level [U] would strengthen the 
bounds from the neutron [H] to ^ < 3 ■ IQ-^^ ^nd S,w,3- IQ-'^S < (M^/3 • 10^ GeV)^. More 
quantitative statements could be made with lattice-QCD calculations of the EFT LECs. 

Additional information comes from the ratio Aid/dd- For 6, rridlMdl is expected to be 
larger than \dd\, whereas for the dimension-six sources we expect md\A4d\ to be of similar 
size or somewhat smaller than \dd\. For 6, Aid is determined by pion exchange, and we can 
again use the link with isospin violation [14] to find Aid — 2.0 ■ 10~^^ e fm^. An upper 
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TABLE I: Orders of magnitude for the deuteron EDM (in units of em^ ), the ratio of deuteron- 
to- neutron EDMs, and the ratio of the deuteron MQM and EDM (in units of tu^^), for sources 
of effective dimension up to six. 



Source 


9 qCEDM qEDM gCEDM 


ma dd/e 

dd/dn 
rridMd/dd 


-j^ QCD -j^ 
^^^QCD 1 1 



bound on Aid can therefore constrain 6 without relying on an estimate of short-range physics 
via the size of the chiral log, which is necessary when using dn [B]- Moreover, if m^lA^fil is 
found to be much smaller than \dd\, the source would likely be qEDM. This shows that a 
measurement of A^^, in addition to dn and dd, would be very valuable, and as a consequence 
its feasibility is beginning to be investigated [25] . 

The deuteron EDM and MQM were calculated previously in Ref. [26]. Since these calcula- 
tions did not use the chiral properties of the fundamental sources, the pion-nucleon 
interactions were assumed to be all of the same size. When the dominant source is the 
qCEDM, their results agree with ours. The advantage of our EFT framework is that it 
has a direct link to QCD by exploiting the chiral properties of the dimension-four and 
-six operators. This is demonstrated by the g2Nii3T^N interaction used in many previous 
calculations, which due to its chiral properties only comes in at higher order for all 
sources [14|. Consequently, for the qCEDM, the ratio of dd to A^^ depends at LO only on 
the ratio Qi/qo, which can be measured independently: cji could be inferred from dd, and 
go in principle from another observable, such as the proton Schiff moment [16J or the ^He 
EDM [27]. In addition, the power-counting scheme allows a perturbative framework with 
analytical results that can be improved systematically. Under the assumption that higher- 
order results are not afflicted by anomalously-large dimensionless factors, the relative error 
of our results should be •j/Mnn ~ 30%, as was explicitly verified for the charge FF [21]. 
Our estimates for dd are consistent with those from QCD sum rules [28j . 

In summary, we have investigated the leading-order, low-energy electric-dipole and 
magnetic-quadrupole form factors of the deuteron that result from the 6 angle, the quark 
electric and chromo-electric dipole moments, and the gluon chromo-electric dipole moment. 
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While for qCEDM we expect \dd\ to be larger than \dn\ by a factor 0{MQQj^/mT^MNN), for 
the other sources we have shown that da is given by the sum of d„ and dp. Furthermore, 
the SM predicts md\M.d\ to be larger than \dd\, whereas beyond-the-SM physics prefers 
iTT-dl-Mdl smaller than, or of similar size as \dd\- EDM and MQM measurements are therefore 
complementary. 
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